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Abstract 



We study an exact swiss-cheese model of the Universe, where inhomoge- 
neous LTB patches are embedded in a flat FLRW background, in order 
to see how observations of distant sources are affected. We find negligible 
integrated effect, suppressed by (L/Rh)^ (where L is the size of one patch, 
and Rh is the Hubble radius), both perturbatively and non-perturbatively. 
We disentangle this effect from the Doppler term (which is much larger 
and has been used recently [1] to try to fit the SN curve without dark 
energy) by making contact with cosmological perturbation theory. 
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1 Introduction 



The fact that our universe seems to be accelerating presently has been one of the most 
remarkable and baffling findings of recent cosmological observations. The origin of 
such a late time acceleration has been a source of intense research (see for example [2]). 
The usual ways of explaining it are either to postulate the presence of a dark energy 
component (an unknown component with negative pressure), or to modify gravity. 
The ACDM model has become the most popular choice being the most economical 
model that is consistent with most cosmological observations, but it is safe to say 
that no completely satisfactory or compelling explanation has yet emerged. There 
are two aspects of the phenomenon which make it very hard to explain: the smallness 
problem, which essentially refers to the fact that the scale of new physics required to 
explain dark energy, ~ meV, is much smaller than expected. Secondly, if the equation 
of state for dark energy is so different from ordinary matter, as observations suggest, 
then one is left wondering about their approximate equality at the present epoch. In 
other words why did the dark energy start to dominate "precisely" today? This is 
also known as the coincidence problem. 

The problem of dark energy is so puzzling that it is perhaps worth taking a step 
back and investigate the assumptions based on which its existence has been inferred. 
The most direct evidence of dark energy comes from looking at distant supernovae, 
or more precisely at the relationship between the luminosity distance and the 
redshift z, upto z ^ 0{1). It turns out that, in order to account for the observations, 
one requires a dominant dark energy component with an equation of state ~ — 1 
[2]. Indirect evidence for dark energy comes from CMB which concludes that our 
universe is flat, VLtot ~ 1, with fi^ ~ 0.3 [2]. In other words, non-relativistic matter 
(dark matter plus baryons) only accounts for around 30 % of the total energy budget. 
The rest, 70% of the energy content, can then be accounted for by dark energy. The 
question is: are these conclusions robust or have we made one assumption too many? 

For instance, an assumption in deducing the equation of state for dark energy from 
the supernova curve is that our universe is homogeneous and isotropic, and that we 
can apply Friedmann-Lemaitre- Robertson- Walker [i.e. homogeneous and isotropic) 
cosmology. Indeed, the FLRW model is a very good approximation in the Early 
Universe, as probed by the homogeneity of the CMB (density contrasts in photons 
and matter are of the order 10~^ and 10~^ respectively). However, at low redshifts, 
the density contrast in matter grows up to values of 0{\) and beyond: small scales 
become nonlinear flrst, and then more and more scales enter this regime. Today, 
scales of order of 60//iMpc (which is 2% of the size of the horizoiJl) have an average 
density contrast of order 1. Structures observed with density contrast larger than 0.1 
(so, already in a mildly nonlinear regime) extend to few hundreds of Mpc (around 
10% of the size of the horizon). Since all our deductions about cosmology are based 
on light paths, and the presence of inhomogeneities affects the background [31I11E1E], 
the luminosity distance and the redshift, one can wonder whether their effects can 
in fact account for the observed Dl{z) curve, without resorting to dark energy. One 

■^In this paper we use h ~ 0.7 
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attractive feature of this possibility is that it has the obvious potential of addressing 
both the smallness and the coincidence problems. Hence, recently, such a possibility 
has received a considerable attention and there is an ongoing debate as to whether the 
effect of inhomogeneities can be substantial. The purpose of this paper is to clarify this 
issue using a (semi-)realistic inhomogeneous cosmological model and performing both 
perturbative and non-perturbative analysis, both analytically (wherever possible) and 
numerically. 

Before proceeding further, we should point out that accounting for the luminosity 
distance-redshift curve is not enough, one also has to address CMB observations. This 
is a much more elaborate task which we leave eventually for future, but we only point 
out here that there are several degeneracies involved in CMB measurements which 
one may be able to exploit to fit the data with non-standard cosmological parameters 
(see for instance [8] for such an attempt). 

Let us start by briefly going over the usual arguments as to why it is believed that 
the effect of inhomogeneities is small and cannot account for dark energy. Firstly, it 
is generally argued that the Newtonian potential $, which characterizes the deviation 
from the FLRW universe, is small even today because it is given approximately by the 
ratio {L/RhY ~ 10"^ at the scale of nonlinearity, L ~ 50//i Mpc. Rh ~ 3000//i Mpc, 
denotes the Hubble radius. Since the corrections to redshift or luminosity distance is 
governed by $, one expects them to be small. Secondly, even though light, as it passes 
through underdense or overdense regions is expected to get preferentially more red or 
blue shifted respectively, these contributions may cancel in the end. This intuition is 
corroborated in first order perturbation theory where it is known that the corrections 
to the redshift, except terms which depend on position (linear Sachs-Wolfe effect) 
and velocity (Doppler term) of sources and observer, precisely vanish. So how can 
inhomogeneities play any significant role? 

To answer this question it is crucial to realize that the magnitude of corrections 
due to inhomogeneities on light paths basically depend on two dimensionless quanti- 
ties, the length scale or period of fluctuations governed by the ratio L/Rh, and the 
amplitude of fluctuations, governed hj 6 = {p — {p))/ (p), where p is the energy density 
and (...) stands for a spatial average. Although the latter can be large, the former for 
realistic inhomogeneities is always small and therefore a perturbative analysis in the 
former is always good. It becomes important then to first identify the power oi L/Rh 
by which the corrections are suppressed. The second part is to determine whether 
the "prefactor" , which depends on 6 and can have non-perturbative corrections, could 
circumvent the suppression to give us something significant. These are the two issues 
we want to address in this paper. For this purpose it is convenient to distinguish 
between two different ways that "inhomogeneous cosmology" may work: 

Scenario I, Local Effect: There could be an effect due to our special position with re- 
spect to inhomogeneities. For instance if we are located in an underdense region, then 
local observations (of redshift and luminosity distance) would be affected due to local 
variations in Hubble expansion. In particular what has been advocated in a long list 
of papers (see [9], [TOl [TTl [121 [131 El [13 [I]) is precisely such a model where the local 
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underdensity can account for a larger (than average) local Hubble parameter thereby 
reconciling the discrepancy between nearby and distant supernovae, without dark en- 
ergy. However, in order to make this work, one requires the local underdense region 
to span at least upto the nearby supernovae i.e. upto 2; ~ 0.1 which corresponds to 
300—400 Mpc/h. Although it has explicitly been demonstrated using exact Lemaitre- 
Tolman-Bondi (LTB) solutions of Einstein's equations that such an inhomogeneity 
can in fact accommodate quite well the supernova curve [1] , the length scale required 
is too large as compared to the observed scale of non-linearity and therefore does not 
provide a completely satisfactory resolutioij^. 

The reason, however, that these models can avoid the arguments discussed above 
is two-fold. Firstly, as we clarify in our paper, the dominant local corrections go as 
V$, which is only suppressed as L/Rh unlike For scales as large as a few hun- 
dred Mpc, the suppression is therefore mild. Secondly, for local observations, since 
light only passes through underdense regions, there is no cancellation of effects and 
it is precisely the Doppler term (which is present already in perturbation theory, and 
which would average to zero for the entire universe) that contributes, enabling one to 
account for the discrepancy between the nearby and distant supernovae [1]. We also 
show how using perturbative analysis one can reproduce the same results, so that 
there is nothing puzzling about this "large" local effect. Physically it amounts to 
saying that, if we live near the center of a large underdense region, all nearby sources 
are attracted towards its boundary: therefore they have a collective radial peculiar 
velocity which mimicks a large local value of H. 

Scenario II, Average Effect: light, while travelling though inhomogeneities, does not 
see the '^average Hubble expansion" but rather picks up systematic corrections to red- 
shift and/or luminosity distance. The hope is that such corrections, which are known 
to be present, could be non- negligible. This is a much more ambitious hope, since one 
wants to have a sizeable effect using only realistic nonlinear inhomogeneities (of order 
50 — 60 Mpc/h) to explain deviation of the Dl{z) curve from the EdS (Einstein-de 
Sitter) universe. A priori there are at least three reasons that may make someone 
hopeful. Firstly, although it is generally believed that $ ^ 1 in our universe, except 
very near black holes, to our knowledge there is no rigorous proof. Secondly, in the 
real universe, since voids occupy a much larger region as compared to structures, light 
preferentially travels much more through voids and this could negate the cancellation 
between voids and structures in quantities such as the redshift and the luminosity 
distance. Thirdly, since 6 ~ 0{1) non-perturbative correction in 6 may enhance the 
corrections pj. We want to investigate whether any of these could give a sizeable 
effect. 

We choose to study a swiss-cheese model containing tightly fitted spherical 
patched. Each of these patches contains a spherical void region and a thin dense 

^We can however point out the striking coincidence that the existence of structures of the same 
size could explain the large angle anomalies in the CMB sky [TH] , and that [TT] reports the observation 
of a 25% lack of galaxies over a (SOO/ZiMpc)"^ region. 

^While completing the present paper, we noticed that a similar strategy has recently been followed 
in [7]. However, we mainly focus on the calculation of the redshift, while [7] focuses on the luminosity 



spherical shell (much like the filamentary structure observed in the large scale dis- 
tribution of galaxies) representing structures. Such a model has its short-comings: 
inside the patch the metric is spherically symmetric, which is not the case in our 
real universe, and spherical symmetry is known to produce surprising cancellations 
absent in real worlqj. Note that this feature suppresses any "backreaction effect" 
on the homogeneous regions: they still evolve exactly as a fiat FLRW Universe. For 
this reason such models do not incorporate any back-reaction effect, but only system- 
atic corrections due to light propagation inside the patches (which is usually called 
Rees-Sciama effect). Secondly, we use an LTB metric which has a critical time when 
structures collapse to infinite density, unlike in the real universe, where the structures 
can virialize. 

Nevertheless, the model is semi-realistic and it is a start if one wants to esti- 
mate corrections to redshift and distances coming from inhomogeneities in our real 
universe. What we find is that although a local Doppler effect [i.e. inside a single 
patch) is only suppressed by L/Rh, any average effect [i.e. which survives even when 
one integrates over a whole patch) is much more suppressed and goes like {L/RhY 
confirming qualitative estimates present in literature [T8l [19]. There are then three 
key questions left to be answered, (i) Do these corrections add coherently which can 
perhaps increase the effect? The answer is yes, and in particular if one is interested 
in looking at corrections at 2; ~ 1, then effectively it goes like {L/RhY (but this is 
still too small to explain the SN curve), (ii) What about the 5 dependent prefactor, 
can it be large? Our estimates indicate that the prefactor is C?(l), and it cannot cir- 
cumvent the [L/RhY like suppression, (iii) What happens if one tries to generalize 
this analysis to perhaps non-spherical cases? Is it possible to get a lesser suppression? 
The answer to this question is well beyond the scope of this paper. We just point 
out that first order perturbation theory rules out an 0{L/Rh) correction, but if one 
could obtain a (L/Rh)^ correction in a single patch, then the "coherence effect" along 
with an enhancement due to large density contrasts might lead to relevant effects. 
We leave further exploration of this possibility for future. 

In section |2l we introduce our inhomogeneous swiss-cheese model of universe and 
discuss qualitatively the expected corrections to redshift and distances. In section [3] 
we discuss in detail the propagation of light, i.e. redshift and luminosity distance, 
for profiles where curvature always remains small (although the density contrast can 
be > 0{1)) along radial trajectories. Next, in section HJ we discuss the connection 
between the LTB metric and perturbation theory about an FLRW Universe in order 
to estimate corrections in non-radial trajectories. In section [5] we discuss profiles 
where curvature is large in some regions (and hence one may expect that $ could 
be larger than what is usually advocated). Finally we conclude and we point out 
possible caveats in our model that still leave hope of finding some non-trivial effect 

distance. 

^The very fact that the two metrics can be glued together is a manifestation of this special 
cancellation: only for spherical symmetry, due to Birkhoff 's theorem, the metric outside the patch is 
insensitive to the details of what happens inside, but it is sensitive only to the total mass, exactly like 
Gauss' theorem in Newtonian physics. Another manifestation of this is, for example, the explosion 
of a supernova: it emits gravity waves only if the explosion is non spherical. 
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with more realistic inhomogeneities. 



2 The Swiss-Cheese Model of the Universe 

In this paper we will consider a swiss cheese type of model to describe our universe, 
where the inhomogeneities reside on spherical "blobs" cut out from a flat homogeneous 
Einstein-de Sitter (EdS) universe. We are going to use LTB metrics to describe each 
of these blobs. As we will see later, LTB metrics can describe a wide range of density 
profiles (as a function of the radial coordinate), allowing us to mimick voids and 
structures. In principle one can consider these blobs to be randomly distributed, but 
since we are interested in looking at inhomogeneities of a particular length scale, it is 
particularly useful to keep in mind a tight packing model, where the structures reside 
at the boundary of the blobs, while the insides are mostly voids. This gives rise to a 
universe with continuous two dimensional structures where voids occupy the bulk of 
the space, somewhat similar to what observations suggest [20]. We will however keep 
our analysis as general as possible with the aim of finding the maximal effect that 
can arise when light propagates through inhomogeneous blobs. 

We stress that there are two different effects: Scenario I. This is present when the 
observer lives inside one of these patches and corresponds to a large local correction 
to the Hubble diagram (see [1] for an analysis of how observations can look like in 
this case) , due to peculiar velocities. To investigate these kinds of effects it is suitable 
to place the observer at the center of a patch. Scenario II. The second effect is a 
net correction that a photon experiences when passing through a blob: in order to 
capture this, we put observers and sources outside the blobs, in the FLRW region, or 
equivalently at two diametrically opposite points on the boundary of the patch. 

To have a picture of our swiss-cheese model it is important to realize that one 
can consistently choose a coordinate system where the boundaries of each of these 
blobs are constant in comoving coordinates. No matter (dust) flows in or out of 
the blobs, and with time the only thing that happens is that the physical radius 
of these blobs expands, exactly obeying the EdS expansion which arise out of the 
junction conditions. This fact enables us to consistently paste these blobs onto the 
homogeneous FLRW metric. To understand this procedure better let us start with a 
brief review of LTB metrics. In the subsequent subsection we will discuss the junction 
conditions. 

2.1 LTB Metrics, a Short Review 

The LTB metric (in units c = 1) can be written as 

ds^ = -df + S\r,t)dr^ + R\r,t){de^ + sinHd^^), (2.1) 
where we employ comoving coordinates (r, 6, cp) and proper time t. 



Einstein's equations with the stress-energy tensor of dust, imply the following 
constraints: 



S\r,t) = (2.2) 

where a dot denotes partial differentiation with respect to t and a prime with respect 
to r. The quantity p{r,t) is the energy density of the matter, and G = l/m^^i is 
Newton's constant. The functions E{r) and M(r) are left arbitrary. The function 
E{r) is related to the spatial curvature, while M(r) basically corresponds to the mass 
inside a sphere of comoving radial coordinate r pT| [H]. If M(r) is an increasing 
function, it can always be chosen to be 

M{r) = —My , (2.5) 

by a suitable reparametrization of the radial coordinate and we will adopt this con- 
vention for the rest of the paper. Here Mq is an arbitrary mass scale. 

One can easily verify that, depending upon the sim of E{r), (12. 3p admits solutions 
analogous to the flat, closed and open universes [21] 1 . but for our purpose it will be 
sufficient to only look at the open case, E{r) > 0: 

GM(r), , , , , 

R = ^^(cosh«-l), (2.6) 

GM(r) , , , , 

t-tb{r) = j^^^^^(sinhu-u), (2.7) 

where tb(r), often known as the "bang-time", is another arbitrary function of r. It 
is interpreted, for cosmological purposes, as a Big-Bang singularity surface at which 
R{r, t) = 0. This is analogous to the scale-factor vanishing at the big bang singularity 
in the homogeneous FLRW models. It is easy to see that the function tb{r) becomes 
negligible at late times, t ^ ti,{r), and therefore we can just ignore it and set th{r) = 
for all r. 

For later convenience we enumerate the simplified equations for our model: 

27rr 

R{r,t) = (2.8) 



7rV2 
3A;(r)3/2 



= Mt = — \— (sinhu - u) , (2.9) 



^The reader should not be confused by the terminology which is usual in the literature: for 
instance a model with E{r) > can easily mimick a flat FLRW model. It is also possible in principle 
to have a mixed model with some regions where E{r) < and other regions where E{r) > 0. 



where we have introduced the conformal time r, the "curvature function" 

E(r) 

and 

M^^. (2.11) 
mpi 



2.2 Junction Conditions and Constraints on Curvature 

Before we proceed further it is useful to choose conventions which simphfy both 
analytic and numerical calculations. We choose Mq in such a way that the "coordinate 
density" Mq coincides with the average (EdS) density po, at the present time to: 



Or, in other words 



M* = p„ = ^. (2.12) 



r3 = toM = ^=. (2.13) 
V Gvr 

Essentially such a convention implies that the coordinate r is approximately the 
physical distance between the comoving coordinate point and the center, today. With 
this convention we also find that the present Hubble scale Hq (and hence the Hubble 
radius Rh) of the average model is simply related to the parameter M: 

R],' = Ho^^^ = ^M. (2.14) 



Next, let us look at the curvature function. There is a procedure on how to match 
the blobs to a homogeneous metric consistently [22] , according to which the curvature 
function has to satisfy 

k'{L)=0. (2.15) 

We denote by L the comoving radius of the boundary where the blob meets the ho- 
mogeneous background. In general a given blob maps to a particular FLRW universe. 
Note that FLRW universes are characterized by a single parameter: it could be the 
value of the "curvature abundance" flk at some given epoch (labeled by the value 
of the Hubble parameter, for instance). One can check that if we use the above 
normalization fl2.12p then 

k{L) = ^Qk, for IQkl < 1, (2.16) 

where Qk is the "curvature abundance" of the homogeneous universe at the same 
epoch. In particular, if we want the curvature to vanish, in order to be consistent 
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with CMB, then we would need fi^ ^ 1, and for simphcity we will only look at 
profiles where it exactly vanishes. Thus our profile has the following property 

k{L) = k'{L) = 0. (2.17) 

There is also an additional constraint on curvature that comes from requiring 
continuity at the origin [23] : 

A;'(0) = 0. (2.18) 

It is now clear that as long as we have profiles, or k{r)^s, which obey (I2.17P and 
fl2.18p . we can embed these spherical inhomogeneous blobs inside EdS consistently. 
In particular we can study how light travels through a series of these blobs. 

2.3 Qualitative Discussion of Corrections 

Having chosen the model, we want to study observational quantities of distant sources: 
the redshift z and the luminosity distance D^. Our aim will be to check whether there 
is any non-perturbative large correction to the FLRW predictions, or if we essentially 
recover results from perturbation theory. Here we discuss what are the possible 
outcomes. 

In order to understand the corrections to various quantities due to inhomogeneities 
it is important to realize that there are two important parameters in the problem: the 
curvature inside a patch (which can be parameterized for example by the maximal 
value fcmax of the function k{r)), and the ratio L/Rh ~ ML. While the former can 
most directly be related to the amplitude of density fluctuations (as we will see, the 
density contrast goes like 6 oc /CmaxT^, for small 6), the latter governs the period or 
length scale of fluctuations. During non- linear structure formation 6 ~ but a 

realistic value of the nonlinear length scales gives at most (ML) ~ 0.01. One can push 
this limit to somewhat larger values but still it remains a small parameter. Therefore 
it plays the most important role in determining the magnitude of correction that one 
expects from inhomogeneities. In particular it is of crucial importance to determine 
the power of ML that appears in the correction to quantities such as 6z (the correction 
to redshift w.r.t. the FLRW result) or SD (the correction to the angular or luminosity 
distance w.r.t. the FLRW result). The next important question is to see whether 
the fcjiax dependent factor in the corrections can ever overcome the suppression due 
to (ML). 

For clarity, we will approach the problem in stages. We will start as a first step 
by considering radial photon trajectories passing through profiles where the curva- 
ture always remains "small" (but 6 can be large). As we will see, "small curvature" 
implies that the Newtonian potential $ ~ 5{L/RhY, and therefore is small as it is 
generally believed. Also, radial trajectories necessarily have to pass through both 
voids and structures, which may have cancelling effects on the corrections. So, we 
are possibly underestimating the corrections, but in this special case one can at least 
find analytical solutions for t{r) (time along the photon trajectory), z{r) and D^lr) 



(and hence also Dl{z)), and therefore precisely determine how the corrections to 
these quantities depend on kja^x and ML. We find that there are local corrections 
(Doppler effect) which go like L/Rh and therefore can be significant. In fact these 
are the kind of effects which are exploited in [T] to explain away dark energy with an 
L ~ 300 — 400Mpc//i, i.e. scenario I. However, overall corrections that are relevant 
for scenario //are much more suppressed and go like {L/RhY, which is too small to 
be significant for supernova cosmology. The /Cmax dependent prefactor cannot over- 
come this suppression. This is also born out by numerical solutions of the equations 
of motion. 

Next we check whether there can be any enhancement of the effect by considering 
non-radial trajectories which do not have to pass through "equal amounts" of void and 
structure. Numerical solutions show that one does not find any significant difference in 
the magnitude of correction, and this can be further corroborated using perturbation 
theory. Thus, at least for low curvature profiles, we are tempted to conclude that 
the average corrections are suppressed by [L/RhY, and there is no significant non- 
perturbative enhancement in /cmax which can remotely overcome such a suppression. 

Next, we consider profiles which have some regions with high curvature. Unfortu- 
nately in this regime one does not have complete analytic control. It is not clear, for 
example, whether it is possible to go to a Newtonian gauge where $ remains small 
everywhere, but this is perhaps also the reason to be hopeful. It turns out that it 
is still possible to estimate the corrections that one can expect in such generic sit- 
uations. Unfortunately, it seems that the corrections look very similar to the small 
curvature case. In other words the Doppler contribution goes as L/Rh, but the over- 
all correction is again suppressed as (L/Rh)^. The only difference, it seems, is that 
the non-perturbative effects in /c^ax can give rise to a larger pre-factor. However, this 
enhancement still falls far short of what could be significant for dark energy. 

Finally, we look at a series of patches to check whether the correction in each patch 
adds coherently or not. This turns out to be the case and thus if we are interested in 
estimating corrections for a total distance that correspond to redshifts of which 
is approximately the same as the horizon scale, the average correction gets a Rh/L 
enhancement. This leads ultimately to an overall effect ~ {L/RhY, but still clearly 
too small. 

In the rest of the paper we provide details of the above discussion. 



3 Small Curvature Regime 

To get an intuitive picture of the swiss-cheese LTB model we will first employ what 
we call "small curvature" approximation [1], which allows us to study everything 
analytically. The strength of the approximation lies in the fact that it can accurately 
describe the dynamics even when Sp/p ^ 1 [I], and thereby allows us to study the 
effect of non-linear structure formation on the different physical quantities. We will 
see in the next section that this is also a regime when perturbation theory is definitely 



valid, i.e. one can compute the Newtonian potential and show that it is small. In this 
section we will focus on only radial trajectories, where a cancellation of effects between 
structures and voids may be at work. So, we are perhaps underestimating the effects 
of inhomogeneities but the hope is that we will have an analytical understanding of 
how the corrections depend on the two parameters, ML and fcmax and whether non- 
perturbative effects in the latter can overcome the suppression due to the former. 



3.1 Small-w Approximation and Density Fluctuations 

In f l2.8p and fl2.9l) , if we only keep next to leading terms in u we find [T] 

R{r, t) ^ ^7VV [1 + i?27'r2fc(r)] , (3.1) 
where 7 and R2 are constants: 

-d (3.2) 

One can check that this approximation is valid [i.e. u is small) when 

u^Uq = -fT>yk{r) < 1 . (3.3) 

We will first consider density profiles where this condition is satisfied for all r 
inside the patches. Now, to compute p(r, t) we have to compute R'. From fl3.ip we 
find 

R\r, t) = ^7^2 [1 + R2i't' (rkir))'] . (3.4) 

From ( 12. 4p we have 

p r, t) = . 7T ■ 3.5 

67r(Mt)2 [1 + R2-f^T^k{r)f [l + R2-f^T^ {rk{r))'] 

We note that the correction term in the second factor of the denominator is actually 
proportional to Uq 

R2-f^k{ry = R2UI < 1 , 

and therefore it can be ignored consistently within our "small-u" approximation. 
Thus we have 

pirA) = = , where Air) = irk)' . (3.6) 

^ ' 67r{Mt)^ [1 + i?27'rM(r)] v ; v ; v ; 

We observe that the FLRW behaviour for the density is given by the prefactor 
multiplying the denominator, while the fluctuations are provided by the presence of 



in 



A{r) in the denominator. It is clear that the density contrast (6) defined in the usual 
manner 

p{r,t) - {p){t) 



is controlled by 

One obtains a simple relation 



r,t)=R2Yr'A{r). (3.8) 



5 ' R.l'r^Ajr) 

1 + e 1 + i?27V2A(r) ■ ^ ^ 

When A{r) is close to its maximum value we have a void, while when it is close to 
its minimum, it signals an overdensity. Note that [1], when e <C 1, 5 ~ e and we are 
in the linear regime: an initial density fiuctuation grows at small times as ~ ts, 
in agreement with the prediction of cosmological perturbation theory. When e is no 
longer small the density contrast grows rapidly (and this result is the same as found 
within the Zeldovich approximation [2l]). 



3.2 Radial Photon Trajectories 

As we will see, we are also able to analytically study the trajectory of a light ray, 
which propagates radially inside one patch. 

The first step is to solve the null geodesic equations. In other words, we have to 
obtain t{r) for a photon trajectory converging at say r = ro at time t = to (from now 
on, the subscript O will denote the observer). The equation for this radial photon 
trajectory is simply given by 

dtjr) ^ ^ R'{r,tir)) 
dr ^yl + 2E{r) ' 

The ± depends on the direction of photon propagation. Now we are only going to be 
interested in placing the observer either at the boundary, r = L (and to be definite say 
at 6 = (j) = 0) when investigating average effects {scenario II), or at the center, r = 0, 
while looking at local effects [scenario i). It is then convenient to define a different 
variable r G [—L, L], such that r = r along the interval r E [0,L], (p = 6 = 0, while 
r = —r for the portion of the path r G [0, L], = vr, ^ = 0. Defining k{—r) = k{r) 
for negative argument^, one then finds that the differential equation that is valid 
through out the path for photons which are travelling along increasing values of r is 
given by 

dt _ R'{7,t{^) ^g^^^^ 



dr ^1 + 2E{?) 

The prime now indicates a differentiation with respect to r. In the rest of the paper 
we are going to drop the tilde from r, but the readers should note that when we talk 
about negative values of r, we are really talking about r. 



^Since the derivative of k{r) vanishes at r = 0, fc(r) remains a smooth function. 



n 



Now, we are interested in solving the above equation in the small curvature regime, 
i.e. keeping only upto linear terms in k{r). One can, in fact, find an analytical solution 
(see appendix which looks like 



r = rir(r) + 5r(r) , 



(3.12) 



where tf is the FLRW trajectory 



Tp{r) = To + -l'^{r - To 



" 2- 



(3.13) 



and 





2 



dr kr 




is the correction coming from inhomogeneities. Here tq is an integration constant 
denoting the conformal time at r = tq and we have defined the dimensionless variable 



for convenience. 

It is now instructive to look at all the different terms and estimate them. The 
first three terms all vanish when the source is at the boundary r = —L although, 
inside the patch, these terms are non-vanishing and dominate over the contributions 
coming from the integrals. The contribution of these first three terms can be combined 
conveniently to give 



As we will see later, this corresponds to the Doppler term in perturbation theory 
and in principle can give rise to large correction to the redshift. In fact it is this 
contribution which has been exploited in the "void models" to explain away the 
dark energy problem. However, this effect only lasts as long as one is inside the 
patch and therefore to account for dark energy one has to have a really large void, 
L ~ 300 Mpc/^, which would at least include the nearby supernovae. We leave for 
future to study how likely such an explanation can be [25] . 

It is easy to see what is the suppression in powers oi L / Ru oi the terms in fl3.14p : it 
sufficient to count the powers of r. The Doppler term 03.161) is therefore proportional 
to (ML), but since it originates as a total derivative it vanishes outside the patch; also 
in the language of first order perturbation theory a Doppler term cannot contribute 
to any average effect. Thus although we have verified this for radial trajectories 
in a spherically symmetric model (and for only a single patch), the robustness of 
perturbative arguments suggests that there are no average corrections proportional 
to {ML) in general. 

Next, let us look at possible contributions of the order {ML^ which one can 
naively suspect to arise from the first integral in 03.141) . However, in the special case of 



r = 



Mr, 



(3.15) 




(3.16) 
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radial trajectories at least, this vanishes when integrating from boundary to boundary 
because the integrand is an odd function of f {k{r) is even by definition). However, one 
may hope that in a more general case (relaxing the spherical symmetry or considering 
non-radial trajectories) such a term might give a non-vanishing contribution. We will 
come back to this point later in section |H 

Thus in the end we find that the only non-zero contribution comes from the last 
term and is proportional to (ML)^. One can in fact come up with an upper bound 
for such a contribution (see appendix \^ in our swiss-cheese model: 

— < ^V6^{MLf ^ l.2{MLf . (3.17) 



3.3 Redshift and Luminosity Distance 



In the previous subsection we have found t{r) along a photon trajectory. In this 
subsection we obtain the redshift z{r) corresponding to a source located at r. The 
differential equation governing this relation is given by [9]: 



dz 
dr 



(1 + z)R' 



(3.18) 



In order to solve ( 13.18P we need to compute R'. In the small-w approximation we 
have 

i?' = ^^^[l + 2i?27V(rA;)']. (3.19) 
9r 

First let us focus on computing "local effects" (which show up when the observer is 
located inside the patch, e.g. in the center). Inside the patch the largest correction 
in fl3.19p is the term linear in f and hence we can ignore ~ in (13.181) (since the 
patch is assumed to be much smaller than the horizon). In this approximation, using 
fl3.16p . one can solve for redshift analytically (see appendix iBl) : 



l + z 



-To 


2 




exp 







AttR 
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(3.20) 



Let us try to understand the corrections fl3.20p . We note that, for an open universe 
k > and f < for our trajectory, so that the correction is always positive. In other 
words one always perceives a greater redshift as compared to FLRW. Why is that? 
In an open LTB patch one has to have an underdense region at the center and an 
overdense region near the boundary. Thus let us consider the simplest profile (the 
argument goes through for more complicated profiles as well), which is to have an 
underdense region at the center, followed by an overdense region which matches to 
the average density region at the boundary. In this case, the quantity rk{r) in the 
exponent generically reaches a maximum somewhere in the middle and then falls to 
zero at the boundary. The maximal deviation from FLRW is therefore attained at 



this maximal value which coincides with the average density {[rk{r)]' = corresponds 
to having no corrections in the density function, see (13. 6p ). Physically: starting from 
the center till the average density point the space describes a void which is being 
stretched more than the homogeneous FLRW. Hence we get an additional redshift. 
After that point the space describes a structure which is contracting, which keeps 
compensating for the extra redshift and finally at the boundary precisely cancels it. 

To see why the "local effect" in redshift is significant, let us look at small redshifts 
{i.e. small f). In this case we find 

z ^ -r^(l + 2R2^^kTl) = ZF{r){l + 2R2^^kT^) = ZF{r){l + 2R2ul) , (3.21) 
9rn 



where ZF{r) is the FLRW value of the redshift. Just as an estimate, taking Uq ~ 1, 
we therefore find a 10% correction. If one allows for profiles where curvature can 
be larger, i.e. Mq > 1, one typically gets an even larger correction. In any case it is 
already clear that locally one can get significant effects from inhomogeneities. 

However the situation is very different when one tries to compute the "net" effect 
of the inhomogeneities on redshift as seen by an observer outside the patch, which is 
relevant for the success of scenario II. We notice that, at the boundary r = —L, the 
redshift in (13.201) coincides with the FLRW result since the curvature term vanishes 
at the boundary. So these terms do not lead to any net or average effect, as was also 
found in the case of 6t in the previous subsection. This can again be traced back to 
the fact that the corrections linear in f come from total derivative Doppler terms. 

To find net effects we have to go to higher order terms (in f). The equation for 
redshift (I3.18P in the small-w approximation gives us 

= h + 2i?27V2(ffc)' - kf^] dr . (3.22) 

This, up to linear terms in k, reads 
27r7 



ln(l + z) 



9 



/diT' f f IxiT' 

—^^2PWjdfMfk)'-jd,- 



Tf 



(3.23) 



We note here that for the purpose of estimating net effects we have to place the 
observer and the source at two diametrically opposite points in the boundary, so that 
now the integral runs from —L to L. This essentially guarantees that the correction 
can only go as (ML)^. The linear terms are absent because they originate as total 
derivatives and k vanishes at the boundary, while the quadratic terms vanish because 
the integrand is odd in f. This is obvious for the last two terms in (I3.23p . The 
first term may look more complicated because we find expressions containing double 
integrals (using [371^ ). However, it is easy to see that all these double integrals can 
be converted to single integrals by suitable integration by parts. As a result a typical 
correction term looks like 



f^'lE J dff'^k{f) 



14 



where the subscript "i?" stands for quantities evaluated at the source (emitter) of 
the photons, and m, n are integers, with m > 0. Again when one is integrating from 
boundary to boundary the integral vanishes if m is odd. Thus the largest correction 
arises when n = and m = 2, giving rise to a correction proportional to (ML)^. 

Actually, one can check that even the (ML)^ term vanishes at linear order in 
k{r), consistent with first order perturbation theory (see next section), but for higher 
orders in k{r) this does not happen. Thus the above analysis is instructive to see 
why we will in general get a (ML)^ correction in the average effects, while we get a 
(ML) correction for local effects. Similar arguments can also be applied to estimate 
corrections in 6D as well, yielding similar results. Thus, to conclude, we expect local 
effects in the Dl{z) curve to be suppressed only linearly in {ML), while average 
corrections receive a large cubic suppression. 



4 Non-Radial Trajectories 8z Perturbation Theory 

The radial trajectories that we discussed in the previous section necessarily go through 
structures as well as voids and this, one may suspect, could provide an additional 
mechanism for cancellation of effects. We now want to test whether considering more 
general trajectories leads to a larger effect in 5z or 5Dai or not. This is relevant 
for our universe because voids today occupy a much larger volume as compared to 
structures and therefore one does not expect the photons to see "equal amounts of 
over and underdense regions"!^. 

In this subsection we consider non-radial trajectories which enter the boundary 
not directed along the radius, but at an angle to it. We then keep track of the redshift 
and luminosity distance as the photon exits through another point at the boundary. 
Since the spherical inhomogeneous patches have an overdensity near the periphery, by 
varying the entry angle we can vary the relative amounts that the photon spends in 
over and underdense regions and thereby break any cancellation effects between over 
and underdense regions. In particular there are trajectories where the photon only 
sees an overdensity and therefore we expect to see an enhancement in the corrections. 
The crucial question is: by how much? Unfortunately, we cannot solve the equations 
analytically as we did in the radial case, but we have to resort to numerics. Numerical 
solutions show no relevant variation between corrections in radial versus non-radial 
case (We show the comparisons in figures ([3]) and (j4])). Thus we conclude that the 
effects must still be suppressed by an (ML)^ factor as before. 

In fact, this result can be confirmed using perturbative estimates. As we will 
see, small curvature also implies that the Newtonian potential is small and therefore 
one can hope perturbation theory to be a useful guide. Moreover, these estimates 

^ Since in open LTB models we can only have structures near the periphery, we can only have 
non-radial trajectories which preferentially cross more of structures and not the other way round. If 
we consider a closed LTB model we could engineer the other case. However, it should not matter, 
since the purpose is really to estimate the correction. 

1.^ 



do not rely on trajectories being radial and is equally applicable for non-radial ones. 
Therefore below we clarify the correspondence between standard perturbation theory 
and LTB models in small-w approximation and also estimate corrections that one can 
generally (radially and non-radially) expect in redshift and distances. 



4.1 Newtonian Potential 

Our aim in this section is to try to cast the metric (12.11) in the small-u limit in a 
perturbative form: 

ds"^ = a\T)[-{l + 2$)rfr2 + (1 - 2<l>)dx^] , (4.24) 

and verify that $ ^ 1 and that it is time independent. We will see that this is at least 
possible in the small-w and small- limit. Rather than trying to find the appropriate 
gauge transformations (which is explained in the appendix [C]), our strategy would be 
to compute for the LTB metric the two gauge-invariant combinations that are well 
known in first order perturbation theory. If these two potentials turn out to be equal, 
time independent and small, then it will be clear that indeed the LTB metric can be 
cast into (K2M . 

In the small-M limit we find that the LTB metric is given by 
ds = jfA-dr +f — ] dr +{ — \ r [l + R,^ r k{r)] dn 



Further, assuming E' <^ 1 we have 



ds^ = ^i -dr' 
M2 



2i?27V2 (^Ffc(r)y - 20i?27VA;(r) 



J. — . — 

dr + r 



Here we have defined the new dimensionless coordinate: 



1 + 2R2^\^k{¥) 
(4.25) 



dn' 
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r = —LrM ^ dT = —Mdr . (4.26) 
9 9 ^ ^ 

Also, the prime now refers to differentiation with respect to r. One can think of (14.251) 
as a truncation of the LTB metric keeping only upto linear terms in k{r). Our aim 
now is to find the gauge invariant variables for the above metric. 

Now, in general, any metric containing only scalar fluctuations can be written as 
ds^ = a^(r){-{l + 2ij)dT^ + 2diUjdTdx' + [(1 - 20)% + DijxWdx^} , (4.27) 
where we have defined the operator 

D,^ = d,dj-U,,V\ (4.28) 



Transforming from the Cartesian to the radial coordinate system we have generically 



-(1 + 2^)dT'^ + 2uj\r)drdT + ^1 - 20 + '^S^ dr^ + (^1 - 20 - r^dft 

(4.29) 



where we have defined 



S^X"--- (4.30) 
r 



It is now easy to read off the different potentials corresponding to the LTB metric: 
= -li?27^ i^T^k + T^7k' - 10i?27^^fc) , (4.31) 

and 

S = 2R2-f\^rk' - 20i?27^^A; , (4.32) 
while the potentials ip, uj vanish, and the background scale factor is given by 

a(r) = ^ . (4.33) 

What is observationally important are of course the gauge invariant combinations 
(that we denote by $ and which in linear order look like 



\i/ =ifj + 

and 



1 












) a 




a 




2 / J 


T 



Ya^Xra)r = l{xrr + Xr^), (4-34) 



^ A > ^V72 f Xt\ , 1 // , 2X'\ Xrttr , . 

$ = + ^Vx [^-^I=^+r\X + — + -7^ • 4.35 

6 a \ 2 / b \ r J 2 a 

$ and \1/ coincide with the potentials ip and in the Newtonian gauge fl4.24p and 
for consistency with first order perturbation theory they should be equal [26]. It is 
a straightforward exercise to compute these quantities using (14. 301) . (14. 311) and (I4.32p . 
Indeed they turn out to be equal and one finds 

_ ^ = -$ = 6i?27^ J dr rk{?) < ^A;^ax(ML)2 , (4.36) 

which is indeed small when the length scales of inhomogeneities that we are consid- 
ering are small compared to the Hubble radius. 

Therefore, one can estimate the Rees-Sciama effect coming from perturbation 
theory and check that it is in agreement with our findings, as we will do in the next 
subsection. 
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4.2 Comparison with Perturbative Rees-Sciama Calculation 

It is well-known that there is an overall effect on the redshift of a photon pass- 
ing through a structure (usually called Rees-Sciama effect [IE, [IS])- This has been 
estimated previously by many authors using different techniques (qualitative ar- 
guments [ISI, Void models [27], swiss-cheese models [IH] , LTB models ([2S1 [23 
EO]), second order calculations [31] )• The aim of this section is to review the logic 
behind this calculation, and compare the second order result with our findings. 

First we present a qualitative order of magnitude estimate, based on the smallness 
of the perturbation of the metric. 

Using a weak- field method (in which one assumes the metric (12. ip . and the grav- 
itational potential $ to be small) it can be shown that corrections to the redshift of 
a photon go like: 



where vq^e is the peculiar velocity respectively of the observer and of the emitter 
in the Newtonian frame. The first two terms are the so-called Sachs- Wolfe effect 
(difference of gravitational potential of emitter and observer). The velocity terms 
are due to the peculiar motion of emitter and observer: they represent the Doppler 
effect. The third term is an integrated term, which is present only if the first order 
gravitational potential evolves with time (Integrated Sachs- Wolfe effect). Now, it is 
well-known that in linear perturbation theory (on scales where we can treat 5 linearly) 
$ is time independent, during matter domination. So, the integrated effect has to be 
zero, and the only effect is due to the position and the motion of the observer (and of 
the source) in the Newtonian coordinates. This means that, even if a photon passes 
through an inhomogeneous blob, this does not matter for its final redshift, since the 
only nonzero effect comes from the starting and the arrival points and not from the 
path between. 

Thus perturbation theory tells us that at first order there are no average effects, 
however there can be local effects! In fact, the Doppler terms in fl4.37p give a contri- 
bution which goes as the radial velocity of the emitter, which in turn goes as ~ V$. 
This gives rise to precisely the kind of terms that we found in our previous section 
while discussing corrections to 5t and 5z along radial trajectories in the small u 
approximation. 

Next, let us look at how "net corrections" can nevertheless arise at the nonlinear 
level in 5. This can be estimated as follows [25] : 



where tc is the typical time scale of change of the gravitational potential and At is 
the travelling time inside a structure. Now, At is roughly equal to the size of the 




(4.37) 




(4.38) 



structure L and tc can be estimated as tc = where Vc is the typical velocity inside 
the structure, leading to 

If the structure is still not virialized a good estimate for Vc is Vc = |V$|. Instead, if 
the structure is virialized a good estimate is given by v'^ ~ In both the cases, we 
may express the result in terms of the density contrast S = Using the Poisson 

equation = AttG p6 we get an estimate in terms of 6 and L. 



For non-virialized structures: 



For a virialized structure: 



5Z ^3/2 (L 



(1 + ^) \Rh 



(4.41) 



We can check now, using a systematic expansion of the metric up to second order, 
if we recover the result fl4.40p . We follow the treatment of [21] to compute the redshift 
of a photon in a Universe with only dust, perturbed to second order. The four- velocity 
of the dust is expanded as: 

U^ = \ + v^^^'' + + . . . (4.42) 

This is subject to the normalization condition If^Uf^ = —1. 

The perturbed spatially flat conformal metric has the metric elements: 

(700 = -(l + 2V^«+V^(') + ...,) , (4.43) 
9o^ = 4^ + \u:f^ + . . . , (4.44) 

= (1 - 20(^' - 5., + ^ + \^ + . . . , (4.45) 

wherJ^ Xj-''^* = and the functions ^/''•''^ 't'^^K and xtj^ represent the r-th order 
perturbation of the metric. 

The first order calculation gives the well-known terms: 

SzW = - + ^«^e. - v^i^'e. - /i(A^). , (4.46) 
where e* is the unit vector in the direction of observation and 



hiTE) 



/•TE 

/ dr 




' TO 


_2dT 



(4.47) 
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Indices are raised and lowered using S^^ and Sij , respectively. 



1Q 



and = + + J^^e' - lx\fe'e^. 

This corresponds to fl4.37p . but it is a more general expression, valid in a// gauges. 
As we just saw, the integrated effect is zero. 

Then, we may study the second order expression for 6z. The general expression is 
very complicated and contains a lot of terms. However, we are interested in studying 
the overall effect from boundary to boundary of a single patch. At the boundary all 
gravitational potentials and all velocities are zero (and even second derivatives of the 
gravitational potential). So the result reduces to 



I boundary to boundary -^2(^e) / 



TE 



TO 



2dT 



(4.48) 



where A^^) = ^i^) + 0(2) + ujfe' - Ixf-e'e^. 



Due to the spherical symmetry of our configuration the tensor traceless part of 
the metric, vanishes. 

Then at second order, in the so-called Poisson gauge (a generalization of the 
Newtonian gauge to second order [HT]) one has: 

VV^)^ = -^r(^'^VV-^'>,, + 2T(,% (4.49) 
where 

V^To = -^((VW-^.;i-^"')- (4.50) 

However one may notice, by simple power counting, that the vector term contains 
one power less of spatial derivatives, which means that it is suppressed with respect 
to the scalar contribution. So, the dominant contribution is: 

S^^lar = -2jdrr Qy^'V, - ^To) , (4.51) 

Let's now try to work out explicitly 5z^'^\ using our spherically symmetric potential 
0(r), given by 14. 361 We get: 

^^L = -2/ drr(^l^--'^^To^ , (4.52) 

where 

r dr r 

To(r) = - j j df [{^'{f)f + 2f^"{f)^'{f)] . (4.53) 

Using (14.361) one finds 

^^fLr = -2jdrr (^6Rh\krr - ^^To) , (4.54) 

90 



where 

r df r 

To(r) = -36i?^74 — df [^{kff + 2f^kk'] . (4.55) 

We notice that the largest non-zero term in (14.541) goes hke k'^{MLY corroborating 
the quahtative estimates for non-viriahzed structures (14.401) . We also notice that the 
cubic dependence in ML agrees with our in the previous section in the small-w limit!"]. 



5 Large Curvature Regime 

In the previous sections we have found that, at least in the small curvature regime, 
there is no significant "average" correction to the luminosity-redshift relation. How- 
ever, in general, LTB profiles can have some regions with large curvature. Thus it is 
important to understand at least how the corrections scale with L in the general case. 
Also, here it is not clear anymore that one can go to a Newtonian gauge where $ <^ 1 
everywhere and so, as suggested in the introduction, perhaps we have the greatest 
chance of uncovering a large correction in such set-ups. Also, in this general case we 
can no longer restrict ourselves to linear terms in fc^ax and therefore any potentially 
large non-perturbative correction in /Cmax could show up here. 

To understand some of the generic features of the model let us fall back to the 
defining equations (12.81) and (12.91) . It is easy to see that (13. 3p can be inverted to 
obtain m as a power-series in t: 

u = Uq{1 + C2M0 + C4MQ + . . .) , (5.56) 

where Uq is given in terms of t via (13. 3p . We had computed the coefficient C2 = —1/60 
to obtain the small-w limit. It is instructive to see how this function behaves for very 
large u: in this regime we can approximate 

sinhM,coshM — * -e". 

2 

From (12. 9p we then get 

lim u — >31nno — ln3. (5.57) 

uo — 

Thus, initially u is proportional to Mq, but later it increases only logarithmically. 
One can also obtain a general expression of i?(r, t) as a power series: 

i?(r, t) = ^-K^V (1 + R2UI + R^ul + ...) = ^7r7 W (l + f{ul)) . (5.58) 

We show the behaviour of the function /(mq) numerically in fig. [TJ The FLRW 
geometry corresponds to just including the first term, and the correction, /, can be 

^^The dependence originates because the first non-zero correction in redshift is obtained from 
second order perturbation tfieory. Tlie correction in tlie proper time, 5t is fiowever linear in k as 
derived in the small- li limit, and one can also check this pcrturbativcly. 
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Figure 1: Here we plot /(ttg), obtained combining (l2SD,(l2J]),(l33]) and ([5381) . 



large if Uq is large (since / is a monotonic function, as it can be seen in fig. [T]) . 
Therefore, at first sight, one might also expect reasonably large corrections. Going to 
the large u limit, where one can compute R approximately, such a reasoning indeed 
seems plausible. By substituting (15.571) in (12. 8p we see that 

suggesting an 0{1) departure from the FLRW universe. In particular it also 
becomes unclear whether there is a suitable gauge transformation which can take 
such a metric to the Newtonian gauge with a small Newtonian potential. Thus it 
becomes very important to check whether indeed one obtains large deviations to 
Dl{z) relation. 

5.1 Estimating Corrections 

What we want to do here is to estimate the magnitude of corrections that one can 
expect for generic profiles. In particular we want to find corrections to the photon 
trajectory governed by equations (13.101) and (I3.18p . In the previous analysis what we 
did was to keep the lowest order terms in both k{r) and f which gave us non-zero 
effects. Here we will consider the full power series expansion in the "amplitude" k{r) 
(since curvature can now become large) but still keep the lowest non-trivial power in 
f. Now, for the purpose of estimation, it turns out that we can ignore the denominator 
(the 1/y/l + E term) in (IS.lOp (see appendixiDlfor details). The evolution equation 
then reads 

^^R'^ ll'r' + i?2n7'V2"(rr )' j . (5.60) 

This equation can be solved in much the same way as it was done in the small- 
u case. The results are also almost identical, the linear terms in r vanish because 
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Figure 2: Here we plot Uq^-^^, obtained combining ([23]) , ([221) , (133]) and i^M 



they originate from total derivative terms, and the quadratic terms vanish because 
the integrands become odd. (Instead of terms hnear in k{r) we have to evaluate 
integrals which involve some powers of k{r), but this clearly does not change any of 
the above arguments.) Thus the largest contribution is again proportional to (ML)^. 
By solving fl5.60p one finds the following correction to Tf (see Appendix |D]) 



^ ^ - [^) E ^2n7^Vo2"2n(2n - 1) / df f'k^ 



(5.61) 



One can also obtain an upper bound for the above correction (see appendix [D|) : 
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< {MLf 
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/ maximum 



(5.62) 



At this point it is important to realize that the function /(mq) is a well defined function 
which does not depend on the specific profiles (see figure [T]) and in fact the function 

0.5 (see figure ([2])). Using this we get the following upper 

\St\ 



u 



2d''f 



has a maximum 



bound for 6t 



< 1.5{MLY 



(5.63) 



not much different from what was obtained in the small-n limit (13.171) . We chose to 
study 6t (even if it not directly observable) to make the discussion simple, but we 
expect similar corrections to redshift and luminosity distance. Numerical solution of 
the equation of motion bear this out and this is what we discuss now. 



5.2 Numerical Results 

In this subsection we discuss the numerical solutions for Einstein equations and for a 
light ray in such a background. 
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Figure 3: Here we plot Am (the difference in magnitudes, m = bLogiQDi, between a given 
model and a reference empty open universe) vs. z (redshift). The red long-dashed line is an 
EdS Universe (only dust), the black dotted line is the AC DM "concordance model" and 
the blue solid line is a model with light going radially through one inhomogeneous patch 
(the observer sitting outside the patch). There is no visible correction to the EdS model at 
z > 0.15 (objects behind the patch), since the correction is very small. We also show the 
density contrast (on the same scale) that the photons see along their trajectory (green thin 
short-dashed line). 

First, we have calculated if there is any sizeable effect with a "large curvature" 
profile in a single patch. We have chosen the following profile: 



which satisfies conditions ( 12.17P and ( 12.18p . and where the amplitude ko is chosen in 
such a way that we have a desired density contrast at the present time to- We have 
then glued this single patch to a flat FLRW background. 

In fig. ([3]) we have run the code which calculates the redshift and the luminosity 
distance of the photon for a radial trajectory in this profile, while in fig. (jlj) we have 
considered a photon propagating non-radially in the same background. The deviation 
due to the Doppler effect is visible, but no net effect is well visible in the plots. 

Nonetheless the overall net effect is non-zero and we have considered its depen- 
dence on L/r^or, by running the program with the profile (15.641) for several values of 
L, keeping ko fixed. We have plotted the net correction in redshift (Sz), as a function 
of L in fig. ([5]) which verifies the cubic dependence derived analytically. 

Finally we have considered a "tight-packing" model with several identical patches, 
and a photon propagating radially through all of them. We plot the result in fig. ([6]). 
As one can see the correction is cumulative. 




(5.64) 
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Figure 4: Here we plot Am (the difference in magnitudes between a given model and 
a reference empty open universe) vs. z (redshift). The red long-dashed line is an EdS 
Universe (only dust), the black dotted line is the AC DM "concordance model" and the 
blue solid line is a model with light going non-radially through one inhomogeneous patch 
(the observer sitting outside the patch). There is no visible correction to the EdS model at 
z > 0.15 (objects behind the patch), since the correction is very small. We also show the 
density contrast (on the same scale) that the photons see along their trajectory (green thin 
short-dashed line). 
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Figure 5: In this plot we show the net correction to the redshift {6z) as a function of 
the size of the patch (L), for an observer sitting outside the patch. The triangles are the 
numerical results, while the solid line shows a cubic dependence {Sz oc L^). As expected the 
points follow a cubic dependence. It is interesting to note that, for a patch with 5 « 0.2 and 
a radius of about 200//i Mpc, there would be a 0(10"^) correction for the redshift along 
the line of sight, which would consequently be visible as a secondary effect in the CMB (see 
also [IS]). 
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Figure 6: In this plot we show the net correction to the redshift {6z) for a photon going 
through multiple adjacent inhomogeneous patches. All the patches are identical (with 
physical radius L = 30 Mpc and density contrast 5 I). Each point in the plot corresponds 
to the correction in redshift with respect to the homogeneous case {5z = z — zflrw), when 
the photon comes out from one patch and enters the next one. 

6 Conclusions 

We have studied photon propagation in a swiss-cheese model in which spherical 
patches are embedded in a flat FLRW Universe, both when the curvature of the 
patches is small and large. In both cases the effects are very small if the observer 
sits outside the patch (integrated effect), while they are large if he sits inside (local 
effect). 

We have checked this in the perturbative regime (small curvature), analytically 
and numerically and we have also made contact with perturbation theory around an 
FLRW background. In this way we have identified the contribution of local terms 
(peculiar position and velocity of source and observer) and integrated effects. As is 
well-known, integrated effects appear only at second order, and this is in agreement 
with the fact that we find very small overall effects if the observer sits outside the 
patch. The nontrivial check that we have performed is to study the same thing in a 
large curvature regime. Both by analytical estimates and numerical solutions we do 
not find significant enhancements due to non-perturbatively large curvature terms. 
We also obtained the dependence of the small net effect on the size of the patch, both 
in the perturbative (small curvature) and non-perturbative (high curvature) regimes, 
and we find a cubic suppression in L/Rhor- 

There might still be two logical possibilities to obtain a larger average effect, that 
we have not explored in this paper. The first is to include somehow the possibility that 
structures virialize: the reason one may expect a difference as compared to collapsing 
structures is because the structures no longer contract after virialization, while the 
void regions keep on expanding. One may therefore expect a photon passing through 

9fi 




both structures and voids to experience a comparatively larger redshift. The second 
possibility would be to break spherical symmetry: because of it, in this paper, we 
could not include backreaction effects of the local patches on the "global" FLRW 
evolution, and therefore this leaves open the possibility that significant backreaction 
effect could appear in a more general case. 

Finally we have clarified the origin of large local corrections in LTB (when the 
observer sits inside a patch), as a Doppler effect due to the peculiar motion of nearby 
sources. This effect has been used ([T] and references therein) to argue for the possibil- 
ity of a larger than average local Hubble expansion and thereby explain the supernova 
data without dark energy, as a mismatch between local observations of the Hubble 
parameter and distant supernovae (outside the local patch). This idea works, but it 
requires the existence of a very large underdense patch (the minimal radius of the 
patch has to extend at least up to z < 0.1, which corresponds to 300 — 400/ h Mpc, 
with an average density contrast of roughly 6 = —0.3 ) and it requires that we sit near 
the centeiF^. While this is in contrast with standard structure formation scenario, we 
can however point out the striking coincidence that the existence of structures of the 
same size has been invoked by [16], in order to explain the large angle anomalies in 
the CMB sky. It is also intriguing the fact that [17] reports the observation of a 25% 
lack of galaxies over a (300//iMpc)^ region, which may signal the existence of such 
structures. 
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A "Average" Corrections to r in Small-ix Approx- 
imation 

In the small-M approximation one has fl3.ip - fl3.4l) 

R = + R2ul) = ^7Vr2(l + R^^^T^k) 

/2' = ^7V[l + i?27V(rA;)'] 

and 

i?' = ^7'-[l + 2i?27V(rfc)'] 
y r 

^^We postpone to future work a study of how close to the center we need to be, to be consistent 
with CMB observations. 
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The evolution equation for r is then given by (13. lip 



dr _ f7'M[l + i?27V(rfe)1 
■ 2k{Mr) 



YM[1 + R2Yr\rky - k{Mr) 



(A.l) 



where we have assumed that the photons are travelhng from negative towards positive 
values of r. The above equation leads us to the iterative expression 



TT 



TT 



^ = ^0 + g7 [r- To) + g7 



-^27^ / dr T^(rky — dr kr"^ 
Jo Jo 



tf + 6t, (A.2) 



where Tq is an integration constant denoting the conformal time at r = ro and St 
denotes the corrections to the FLRW trajectory coming from the integrals. Now, 
since we are only interested upto terms linear in /c, we can replace r by Tp inside the 
integrals, so that we have the approximate result 



TT 



St = T — Tp = —7^ 

9 ' 

where we have defined 



and 

h = i?27'^o / dr 



R.,^JdrTUrky-Jdrkr^ 
dr kf'^ , 



_ ^ 2/r 

= 97 ih 



(A.3) 



7r7^r 



I, = 



(fky = R2Yfo / dr 



1 + 2 



7r7^r 



rk + 2—^kr'^ 
9fo 



7r7 
9^ 



kr^ 



, 7^7 



9ro 
df kf — 2 



TT'-f 

mi 



TT'y r 

2N 2 



'-\ 2 



(fky 



dr kr 



In the last step we have used integration by parts and assumed that the observer is 
located either at the center where r = 0, or at the boundary, r = L where k = 0. 
Combining Ji with I2 we have 



TT 



= -7^ I R2l'n 



2^2 




fk + 2— A;f^ + 
Qrn 



7r7 
9^ 



2\ 2 



kf^ - 2 



7r7^ 
9^ 



dr kr 



dr kr 



(A.4) 

The first three terms correspond to local effects which vanish at the boundary. The 
fourth term (integral) vanishes when evaluated from —L to L. So the net correction 
comes from the last term and is indeed very small. In fact one can easily find an 
upper bound as follows: 



\St\ 



2tt 
15 



7 



dr f'k < —^^k^^{MLf . 



(A.5) 



attained by k{r) inside the profile. Now, for small-w we have the obvious bound 



u', = Yr^Hr) < 1 



7 ^max ^ 



^2 ■ 



(A.6) 
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Thus we get an upper bound 
\St\ 27r7^ 



To IStq 

where we have used (12.131) and (13. 2p . 



47J- _ 

df kf'^ < — \/67r(ML)2 ^ 



l.2{MLf 



(A.7) 



B Local Redshift Corrections 



In this appendix we try and obtain the redshift z{r) corresponding to a source lo- 
cated at r inside the patch and the observer at the center. The differential equation 
governing this relation is given by 



dz 
dr 



(1 + z)R' 
VI + 2E 



(B.l) 



To estimate corrections inside the patch, which go ~ r, we can ignore E{r) ~ r^. We 
then have 

^ = -R'dr = -—7' \r-' + 2R2i^T{kr)'] df 
1 + z 9 



27r 2 



TT 



( — ) fk + 2R2Yr{kfy 
9 \tf 



dr . 



where we have used (13.161) . The first term can be integrated in a straightforward 
manner to yield the FLRW result. The second term clearly is suppressed with respect 
to the third term by an extra factor of f and therefore can be ignored. For the third 
term one can use the trick of first integrating with respect to r assuming that r is 
a constant and then putting back the r-dependence of r {i.e. we keep only leading 
powers of r). After all these integrations we find 



ln(l + ^) -lnC= -2 In r^(r) 



4i?27r 4 



9 



7 Tpkr , 



fB.2) 



where C is an integration constant which can be fixed by demanding that at r = 0, 
z = 0. This gives us 

lnC = 21nro, (B.3) 



so that we have 



1 + z{r) 



exp 



-7 Tkr 



9 



(B.4) 
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C From Newtonian to Synchronous Gauge 



Our aim in this section is to try to cast the metric (12. ip in the small-M and small- ii^ 
limit (or equivalently keeping upto linear terms in k) 

ds^ = {-dr^ + [1 + 2i?27V2 (FA;(?))' - 20i?27'^A;(r)] d^^ + [l + 2R2-f^T^k{r)] dQ^} 

(C.l) 

in a perturbative form: 

ds^ = a\T)[-{l + 2$)rfr2 + (1 - 2<!>)dx^] , (C.2) 
where $ ^ 1. We will see that this is possible by a suitable gauge transformation. 
In general, any metric containing only scalar fluctuations can be written as 
ds"^ = a2(r){-(l + 2ip)dT^ + 2diUjdTdx' + [(1 - 2<p)5,j + Dijx]dx'dx^} , (C.3) 
where we have defined the operator 

D,j=d,dj-U,,V\ (C.4) 

The linear gauge transformations which leave the form of the metric and the scale 
factor invariant is given by 

i, = i^-Cr--C (C.5) 

a 

U = U + C + Pr (C.6) 

= <^-^V^/5 + ^C (C.7) 

X = X + 2P, (C.8) 

where the tilded variables are the new coordinate system, (3 and ( are two scalar 
parameters that define the gauge transformation, and the subscript r indicates partial 
differentiation with respect to the conformal time r. 

The strategy that we will adopt is to find appropriate gauge transformations which 
take us from the Newtonian gauge l\C.2\} to the synchronous (LTB) gauge. First of 
all, in order to deal with the LTB metric we have to work in the radial coordinate 
system. To convert the different terms in (IC.SP in the radial system we have to use 
the identities 

= ^ , (C.9) 
r 

and the differential relations 

Xidx^ = rdr , and dx^dxi = dr^ + r'^dVt^ , (C.IO) 

Then we find that when uo , are only functions of the radial coordinate, the terms 
involving them in the metric ( 1C.3I) are given by 

2diuj{r)dx' = 2uj'{r)—dx' = 2uj'{r)dr , (C.ll) 
r 



and 



;i - 2(t))6ijdx'dx^ = (1 - 2(j)){dr' + r'dO.') 



(C.12) 



while the term involving ip remains the same. 



To compute the last term in flC.3l) involving x we need to evaluate the operator 
Dij. Straight forwardly we find 



didjx{r) 







+ 6ij — 










r 



and 

so that 

Dijxdx^dx-' 



r 



r 



1 

3' 



where we have defined 



r 



(C.13) 



(C.14) 



(C.15) 



(C.16) 



Thus putting everything together we have 



ds = a (t) 



1 + 2'4})dT^ + 2u\r)drdT + ( 1 + 20 + ) dr^ + ( 1 + 20 - ) r'^dVt^ 



Now, for the Newtonian gauge, we have 

■0 = = $ and X = = . 



(C.17) 
(C.18) 



In order to go to the synchronous gauge we demand ip and uo in (IC.SP to vanish, which 
gives us 

c, + ^?ic-* = o ^ c = 5^, (C.19) 

a 3 
(here we ignore a possible integration constant), and 



C + /9r = 



/3 



6 



(C.20) 



where the last term arises as an integration constant. Under these gauge transforma- 
tion the non-zero potentials become 



X = + 2/3 
and 



$-2 



|^ + 2/3o(r) ^ 8 = -^- 



1 <I)'t2 9 R' 1 
?> ' T 2, 18 9 r 3 r 3^° 



+ 2(/3^'-^), (C.21 



(C.22) 



We now have to check whether by appropriately choosing $(r) and /3o('") "we can 
obtain S{r) and corresponding to the LTB metric. By comparing the LTB metric 
(IC.ll) with the metric in the radial coordinate system (IC.ITP and using (1C.22P we find 
that this amounts to finding a {$(r), /3o(r)} such that the following two equations 
are satisfied: 



-S 



— $ + 2/3^' - -r2$" = 2i?27V' {rk{r})' - 20i?27'?^A;(?) , (C.23) 
3 3 



and 
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y 



(3' 1 r2$' 

$ + 2^ — 

r 3 r 



The latter equation implies 



3 r 



/9o = 3*^: 



and 



(C.24) 



(C.25) 



(C.26) 



One can now easily check that the same solution set {^{r) , /3o(r)} also satisfies 

(ESI. 



D Estimating Corrections to r in the General Case 



The purpose in this section is to try and generalize the analysis done in the small 
curvature regime so that we can obtain an estimate for the correction for general 
density profiles. The general expressions for R and R' can be written as a power 
series in Uq, or k{r): 



and 



TT 



R' = -7^2 
3 ' 



l + X]^2n7'V"(rfc")' 



(D.l) 



The evolution equation for r is then given by fl3.1ip and fID.ip 
dr f7Ml + ErWv2"(fA;")'] 



TT 



r7 



dr y/l + 2kF^ 9 

where the coefficients Sn are defined by the expansion 

oo 



1 + ^2n7'"r2"(rA;")' + ^ 



^2n 



[D.2) 
(D.3) 



From (1D.2P we find the iterative expression 



T 



Tf 



6t 



Tf 



~9" 



1 ^ 1 ^ 



[D.4) 



Now except for tlie n = 1 term, all other ones in the {Sn} series in flD.2l) are 
suppressed by a factor or more and therefore we can ignore them; as we have 
explained before f <C 1 and higher powers of f will be suppressed. Further, the 
n = 1 term has already been discussed in the small-w case (finding that it goes like 
(L/RhY), so we do not need to keep track of it. Thus we are left evaluating the {Rn} 
series in (lD.2p . To get the leading non-zero correction, it is sufficient to keep only 
upto 0(f^) terms in the expansion of r^" that needs to be substituted in (lD.2p . This 
is straightforwardly obtained from the zeroth order (FLRW) result: 



T ^To 



9 



r 



2n ^ _2n 
~ Tn 



1 : h n{2n — 1] 



9ro 



9ro 



where we have neglected once again higher powers of r. So, St now reads 



2 oo „ 



1 : h n[2n — 1 



9ro 



9ro 



irk 



n\l 



(D.5) 



[D.6) 



We are only interested in estimating the correction when the photon goes from 
boundary to boundary. As in the small-w case, the first term, being an integral of a 
total derivative, vanishes when the integral is evaluated from a boundary to another 
boundary. Next let us try to evaluate the contribution from the second term. It is of 
the form 



rM{rk'')'dr = M 



2un 



r k 



dr rV' 



The first term again vanishes at the boundary or at the center. Moreover, since the 
integrand in the second term is odd, the integral also vanishes when evaluated from 
boundary to boundary. 

Thus we can only hope to obtain a non-zero contribution from the third term in 
f lD.6p . Let us therefore look at the integrals 



f^ifk'^ydf 



f^k'' -2 df f^k"" 



[D.7) 



Now, the integrand in the r.h.s. is even and therefore it does not vanish. This is the 
first non-zero correction coming from the inhomogeneities and gives us 



6t 

To 



(2 \ 3 CO „^ 
^ J J2 W"ro'"2n(2n - ^) J ^df f'k^ . (D. 



One can try to find an upper bound for the above expression as follows: 

J dff'k-<^k:^{MLf. (D.9) 



Thus we have 



To 



In order to estimate the sum, first we observe that the sum can be written in 
terms of a function f{uo): 

1 "Mq 



Therefore we have 



(D.ll) 



^ < {MLf (u''^ 
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